Superconducting Magnetization above the Irreversibility Line in Tl2Ba2Cu06+<5 
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Piezolever torque magnetometry has been used to measure the magnetization of superconducting 
Tl2Ba2CuOe+a. Three crystals with different levels of oxygen overdoping were investigated in mag- 
netic fields up to 10 Tesla. In all cases, the magnetization above the irreversibility line was found to 
depart from the behaviour M ~ \n(rjH C 2/H) of a simple London-like vortex liquid. In particular, for 
a strongly overdoped (T c — 15K) crystal, the remnant superconducting order above the irreversibil- 
ity line is characterized by a linear diamagnetic response (M ~ H) that persists well above T c and 
also up to the highest field employed. 

PACS numbers: 74.25.Ha, 74.60. Ec, 74.72.Fq, 07.55. Jg 



I. INTRODUCTION 

A distinguishing feature of the high-T c cuprate super- 
conductors is that a finite resistance (in the zero current 
limit) appears at a magnetic field well below that required 
to restore the full normal state resistance. Experimen- 
tally, this resistance onset field is found to coincide with 
the so-called irreversibility field H ar , above which mag- 
netic irreversibility vanishes. A broad superconducting 
to normal transition is thought to reflect the existence of 
a mobile vortex liquid between H m and H C 2 .a Although 
diffraction evidence for a genuine vortex liquid has yet to 
be reported, the term will be used here to designate the 
state of the vortex assembly between these two character- 
istic fields. Physically, the existence of a liquid is thought 
to be a consequence of large anisotropy and a short su- 
perconducting coherence length, both factors weakening 
the vortex-vortex interaction. In one of the many avail- 
able scenarios,El this weak flux lattice is thought to be 
"melted" at H lrr by strong thermal fluctuations. 

An open question is how this behaviour evolves 
as T c is lowered and the coherence length increases. 
T^b^CuOg+a is an attractive compound for an investi- 
gation of this issue: with an optimal critical temperature 
of around 90K, it exists in a comparatively simple single 
layer structure, and strongly overdoped tetragonal single 
crystals can be produced by introducing relatively little 
excess oxygen (T c — > for 5 = 0.1 )aU In some respects, 
the overdoping that is possible in Tl2Ba2Cu06+<5 is 
unique: for example, overdoping the L^-zSr^CuC^ sys- 



tem leads to high levels of disorder. In Tl 2 Ba 2 Cu06+5, 
the disorder associated with doping resides in or between 
the TIO bilayers, so disorder-related carrier scattering 
in the Cu02 planes is relatively weak; in-plane resistiv- 
ities just above T c can be less than 10 /lilcm, even in 
strongly overdoped material.u A variety of experiments 
have been performed on tetragonal single crystals and 
epitaxial thin films. For example, the superconducting 
gap symmetry has been assessed in tri- and auadri-crystal 
thin filmscrQ and in microwave experiments^ the low-field 
vortex properties-pf the vortex lattice hane been studied 
by magnetizatioro and Bitter patterning ,E3 and a number 
of normal state properties have also been investigatedBEJ 

In Tl2Ba2CuC>6+<5 crystals with high values of T c , re- 
sistance measurementsEj have shown that the normal- 
superconducting transition is broadened, as observed in 
other cuprate materials. However, for strongly overdoped 
material with-Jow T c values, sharp resistive transitions 
are observed.cl These are similar in appearance to the 
transitions encountered in conventional superconductors, 
where it is known that Hi„ ~ H C 2a This suggests that 
lowering T c by overdoping somehow eliminates the flux 
liquid from the vortex phase diagram. 

On the other hand, the temperature dependence of the 
apparent upper critical field deduced from the overdoped 
resistivity data is strikingly different from that expected 
in a conventional superconductor: it exhibits an upward 
curvature from T c = 15K down to the lowest tempera- 
tures studied (12mK) and shows no sign of saturation 
for T — ► 0.B Qualitatively similar observations have been 
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reported on thin films of I^S^CuOeEl and other hole 
doped cuprates. This anomalous apparent H C 2 in the 
hole doped materials has attracted substantial attention. 
Theoretical proposals fall into two broad classes. One of 
theseli-a'lLa regards the existence of residual superconduct- 
ing order above the apparent H C 2 as an essentialj-aspect 
of the phenomenon, while the other approachE3~E3 con- 
siders the transition seen in the resistivity as an accurate 
estimate of the location of the transition from the super- 
conducting to the normal state. 

The sharpness of the transition and the absence of a 
measurable magnetoresistance at higher fields both argue 
against the existence of residual superconducting order. 
However, it is difficult to estimate quantitatively the con- 
tribution of such order to the conductivity, so transport 
measurements by themselves cannot decide this issue. 
Studies of the specific heat or magnetization would be 
more definitive, but single crystals of high quality tend to 
be small (mass < 10~ 5 grams) and the experiments ccfc 
respondingly difficult. Nonetheless, Carrington et al.E3 
have recently succeeded in measuring the specific heat 
of a small single crystal with T c — 17K. Although an 
applied field dramatically reduces the size of the super- 
conducting specific heat anomaly, the results imply that 
remnant superconducting order persists well above the 
apparent H C 2- 

In summary, the information available on the vortex 
liquid in T^I^CuOg+a presents a seemingly conflict- 
ing picture. As the transition temperature is reduced by 
overdoping, specific heat measurements suggest that a 
wide vortex liquid region persists. By contrast, the sharp 
resistive transitions suggest that the liquid region be- 
comes very narrow in highly overdoped material. Magne- 
tization offers the most direct probe of vortex behaviour, 
but previous SQUID efforts to measure this quantity on 
highly overdoped piaterial have been hampered by inad- 
equate sensitivity.EJ 
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FIG. 1. Predicted angular dependence of the torque (in ar- 
bitary units) according to Eq. (|l|), with H/rjH C 2 =0.1. Ori- 
entations of the applied field H , the magnetization M, and 
the c-axis are shown in the inset. 



We have employed a recently introduced variant of 
torque magnetometry - using piezoresistive microcan- 
tilevers - to reach a working resolution of 10 _13 Am 2 in a 
field of 10 Tesla. This is about two orders of magnitude 
more sensitive than is available from standard SQUID 
techniques, allowing us to explore the magnetization of 
Tl2Ba2Cu06+5 over a wide doping range. We find that 
a diamagnetic torque persists above H m for all the crys- 
tals examined. However, the underlying magnetization 
departs drastically from that of a London liquid. 



II. TORQUE PRINCIPLES 

For a material in which the magnetization M is uni- 
form, the magnetic torque density (torque per unit vol- 
ume) is given by r = M x H, where H is the applied 
magnetic field. If H is applied along a symmetry axis 
of a single crystal, the equilibrium magnetization is par- 
allel to H, and the torque is zero. The magnetic re- 
sponse of tetragonal high-T c materials to off-axis fields 
is largely controlled by their superconducting anisotropy, 
7 = (nic/nia) 1 / 2 . Here, m c and m a are the Ginzburg- 
Landau superconducting effective masses for pair motion 
along the c-direction, and in the Cu02 planes, respec- 
tively. Both the magnitude and direction of M depend 
on 7 and on the angle 9 that the field makes with the 
c-axis. The anisotropy of all the crystals studied in this 
work was large (7 > 50), simplifying the interpretation 
of the data. 

The torque on an assembly of vortices in a high- 
Tc superconductor exhibits a characteristic angular 
dependence.^ Its origin may be understood by referring 
to Fig. [|. For large 7, the magnetization lies very close to 
the c-axis, as shown in the inset, while its magnitude M 
depends only on the component of the applied field along 
that axis. The supposition that M is fully determined 
by the effectives-field HcosO is a fundamental result of 
scaling analysiscB in the large-7 limit of the anisotropic 
Ginzburg-Landau model. 

Under these circumstances, and 

for H c i <C H |cos6>| <C H C 2, the equilibrium vortex mag- 
netization can usually be well approximated by the Lon- 
don result M = (<J>o/87i7i A 2 ) In (r)H C 2/H |cos0|), where 
A and H C 2 are the in-plane penetration depth and the 
upper critical field along the c-axis, respectively, and 77 is 
a numerical parameter of order unity.E3 For < 9 < ir/2, 
the magnetic moment lies along the negative c-axis, and 
the magnitude of the vortex torque density is given by 



Tv (6) = MHsint 



§ H sin 6 



In 



VH C 



H cos( 



(1) 



For 9 > 7r/2, the magnetization points along the posi- 
tive c-axis, and the sign of the torque is reversed. At 
9 = 7r/2, the magnetic field crosses the CUO2 planes, and 
the superconducting screening currents abruptly change 
direction. 
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A more complete description of the vortex torqueO 
indicates that Eq. is correct to ±2% for 7 = 50 and 
for field orientations lying more than 5° from the CuC>2 
planes. A typical predicted torque characteristic is shown 
in Fig. |l|. The drastic increase of r at high angles is 
due to the logarithmic term in Eq. (|l|): as the effective 
applied field decreases, the diamagnetic magnetization of 
the crystal increases, reaching a sharp peak as the crystal 
enters the Meissner state at H c \. 

In traditional magnetization measurements, the irre- 
versibility field is probed by cycling the magnetic field up 
and down and recording the value at which irreversibil- 
ity develops. Since the irreversibility is anisotropic, the 
symbol H m is reserved here to denote the c-axis irre- 
versibility field. For reasons discussed in the Appendix, 
-ffhr was measured here in a different manner: the field 
was fixed at some value H , while the angle 0, and with 
it the effective field Hcos9, was cycled. Angular hys- 
teresis is observed above some angle 6 1TV . Assuming a 
large anisotropy, H lrr can then be deduced from directly 
measured quantities, viz. 



H cos 9i T 



(2) 



piezolever was employed to compensate background sig- 
nals, using a Wheatstone bridge circuit driven by a float- 
ing 77 Hz AC current source. The current (50/iA through 
each lever) increased the sample temperature by less than 
0.1K, at all temperatures. 

The two levers were mounted closely together on a ro- 
tation stage inside a pumped 4 He cryomagnetic system 
capable of reaching 1.3K in a 15T field. The magnet was 
operated in persistent mode, and its field inside the rota- 
tion stage was homogeneous to ~ 0.1%. The cantilevers 
were either in 4 He exchange gas or immersed directly in 
liquid helium. 

Both calibration and uncertainty estimation are more 
involved than for traditional torque methods.E3 The main 
source of random uncertainty comes from variations in 
temperature. A systematic background signal originates 
from the varying gravitational torque on the crystal when 
the sample stage is rotated. However, the largest mea- 
surement uncertainty stems from the intrinsic magnetore- 
sistance of the levers. Over most of our measurement 
range, the torque density data presented in the next sec- 
tion are estimated to be reliable to ±10 N/m 2 . For a more 
detailed discussion of the measurement uncertainties, the 
reader is referred to Appendix [Br 



III. EXPERIMENTAL 



IV. RESULTS 



The critical temperatures of the three Tl2Ba2CuOg+5 
crystals used in this work were established as the tem- 
peratures above which the low-field superconducting vor- 
tex torque signature disappeared. The values obtained, 
T c = 15K, 25K, and 85K, agreed with previous estimates 
from resistivity measurements for the T c = 15K sampletl 
and with thc-spccific heat anomalies observed in the other 
two crystals! 2 !! within experimental uncertainty (domi- 
nated by the superconducting transition widths) of a few 
percent. 

All three samples were grown using a self-flux method 
in alumina crucibles.cS As-grown crystals are overdoped, 
so the desired doping level was established by low tempe*-. 
ature (T < 400°C) annealing in various atmospheres .Ea 
For our torque experiments, the samples were then cut 
into platelets with typical dimensions 160 x 80 x 10 /im 3 . 

Torque measurements were performed using piezore- 
sistive microcantilevers. Originally designed for atomic 
force microscopy,ca their use as sensitivej-tprque sensors 
was pioneered by_Elossel and coworkers. Ej Commercial 
silicon piezoleversEll were employed in this work, their di- 
mensions being 170 x 50 x 5 fim 3 . A boron-doped path 
implanted on the upper lever surface has a resistance of 
about 3kf2 at 4K. On applying a torque, the resistance 
changes by a small amount AR. A typical number for 
the response is AR/R ~ 10 7 per Nm. 

Using an XYZ micropositioner, the crystal was epox- 
ied to the end of a cantilever, its Cu02 planes coincid- 
ing with the flat surface of the lever. A second (empty) 



A. Anisotropy and Irreversibility 

In principle, the superconducting anisotropy 7 may be 
determined from the angular dependence of the torque 
when the field lies close to the Cu02 planes. In prac- 
tice, the finite mosaic spread of the crystal broadens 
the high-angle characteristic. For this, and other rea- 
sons, the data can only be used to establish a lower 
bound for 7. .-These difficulties have been discussed in 
the literature.!^ 2 ] Nonetheless, for all three crystals, we 
estimate that 7 > 50, justifying the high-anisotropy ap- 
proximation described in section O. 

Fig. U shows the angular dependence of the torque den- 
sity observed for the T c — 85K crystal at T — 35. 3K in a 
field of 1 Tesla. (Note that in this and subsequent fig- 
ures showing angular dependences, the data have been 
taken over two angular quadra nts and have then been 
symmetrized according to Eq. (B2), as discussed in the 
Appendix). The irreversibility angle m , identified on 
this Figure, allows the associated irreversibility field H nl 
to be estimated using Eq. (|J). iJj rr (T) was obtained in 
this manner for all three crystals and is shown in Fig. |^. 
The data for the T c = 15K crystal are in good agreement 
with previous reports of the "foot" of the resistive tran- 
sition in a magnetic field£l Also, the irreversibility fields 
obtained with the angular sweep technique agreed with 
control runs at selected temperatures that used a simple 
field sweep with H pointing close to (but not along) the 
c-axis. 
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FIG. 2. Angular dependence of the torque density for the 
T c = 85K crystal in an applied field of 1 Tesla and at a tem- 
perature of 35. 3K. The direction of angular sweep is indicated 
by the arrows. Irreversibility occurs above the angle identified 
as #irr- The complex behaviour inside the irreversible region 
will not be discussed here. 



Note that the irreversibility lines for the T c = 15K 
and the T c — 85K samples almost coincide when plot- 
ted against reduced temperature. We believe that this 
is an intrinsic property of the material - it should again 
be stressed that the samples used in our study are of 
high quality as indicated by their sharp specific heat 
anomalies. Ell Also, the lack of an unambiguous trend of 
the irreversibility fields with T c is again in accordance 
with resistivity measurements across large parts of the 
doping rangeJlj Still, this behaviour is quite surprising 
and warrants further investigation, as one might naively 
rather expect a monotonic change with doping. 
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FIG. 3. The irreversibility field along the c-axis, Hi rT , as 
a function of reduced temperature, T/T c , for the T c = 15K 
(open triangles), 25K (open circles), and 85K crystal (full 
circles) . 



B. Temperatures below T c 



Torque density data for all three crystals, at 
T/T c = 0.85 and H = 5T, are shown in Fig. |. At a 
reduced temperature of 0.85, a field of 5T along the c- 
axis is much greater than H m for all three crystals (see 
Fig. ^ . Irreversibility is therefore confined to angles very 
close to 90°. 

There is a sharp quantitative distinction between the 
three cases. The torque densities for the T c = 15K, 
25K, and 85K crystals stand in the approximate ratios 
1 : 10 : 100. The angular dependences for the T c = 25K 
and 85K crystals are in at least qualitative accord with 
Eq. (|l|) (see Fig. |l]). However, the angular dependence 
for the T c — 15K crystal is quite different, containing a 
strong contribution that appears to vary as ~ sin 29. 
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FIG. 4. Typical angular dependence of the torque den- 
sity for the crystals investigated in this work. In all cases, 
T/T c i± 0.85 and H = 5T. The solid curves are the fits ac- 
cording to Eq. . 
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We therefore tried fitting all the data to an expression 
of the form 

t(6) = Ax sin 29 + A 2 sin 9 + A 3 sin 9 In |cos 6\ (3) 

for < 9 < 7r/2. For 9 lying in the second quad- 
rant, the torque is given by t(9) = —t(tx — 9), in con- 
formity with its antisymmetry around the afr-plane. 
The sum of the A 2 - and ^3-terms is equivalent to 
— A3 sin 6* In (r]H C 2/H |cos#|), with an apparent critical 
field T]H C 2 = -ff exp(— A2/A3), and therefore corresponds 
to the pure vortex torque expected from Eq. ([!]). The 
linear decomposition in Eq. (|) provides the most con- 
venient representation of our data in that it facilitates 
error estimates and in that the fitting coefficients have 
the same units (N/m 2 ) as the torque density itself. Ai, 
A2, and A3 were adjusted to give a least squares fit to 
the data at each field and temperature. A central re- 
sult of our work is that Eq. (^) provides an excellent 
empirical fit to all the data (see, e.g., Fig. ||). The fit 
parameters for all temperatures, fields, and samples ex- 
amined in this study have been tabulated in Table [| and 
therefore provide a comprehensive representation of our 
experimental results. In all cases, the fit is within experi- 
mental uncertainties in the region < 9 < 9- m , where the 
irreversibility angle #i rr can be found using Eq. (Q) and 
the f/irr data shown in Fig. 0. 
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FIG. 5. Temperature dependence of the normalized sin 29 
coefficient A1/H 2 = -Ax/2/io, for the T c = 15K and 25K 
crystals, at H = 3T (circles), 5T (diamonds), and 10T (tri- 
angles). (The H = IT data has been omitted as it is af- 
fected by the onset of irreversibility close to the afe-plane be- 
low T c .) The magnetoresistive error is estimated to be less 
than 0.5N/m 2 T 2 at all temperatures, cf. the Appendix. 



C. Temperatures above T c 

For all three crystals, both the coefficients A2 and A3 
in Table [l|, corresponding to the superconducting vortex 
torque, extrapolate to zero at T c , but this is not the case 
for the coefficient A\\ for the T c = 15K and 25K crystals, 
the sin 29 term could be extracted for temperatures well 
above T c . However, for the T c — 85K crystal, the contri- 
bution of the sin 29 term is very small compared to the 
magnetoresistive background (see Appendix |b|), and no 
systematic data could be obtained above T c . 

The sinusoidal torque A\ sin 29 in Eq. (||) corresponds 
to a magnetization M that is proportional to the effective 
field H cos 9 along the c-axis, with a linear susceptibility 
of x = — 2/ioAi/iJ 2 . (More precisely, if the additional 
effects of small normal state paramagnetism are taken 
into account, this quantity is the susceptibility anisotropy 
Ax — Xc — Xab ) For Eq. to consistently describe the 
torque data at different applied fields, the coefficient A\ 
therefore has to be proportional to H 2 . 

Fig. H shows data for the temperature dependence of 
Ax/H 2 for both the T c = 15K and 25K crystals at differ- 
ent fields and is indeed in basic agreement with A\ ~ H 2 . 
Deviations from an ff 2 -dependence evidently occur at 
temperatures below T c , probably due to slight systematic 
variations between the runs at different applied fields, 
with a comparatively large effect on the outcome of the 
non-orthogonal three-parameter fit in Eq. (^|) . 



V. DISCUSSION 

Data for the angular dependence of the torque 
in T^B^CuOg+s over a wide range of tempera- 
ture, field, and doping are parametrized in Table Q. 
We now discuss the magnetization M, and its de- 
pendence on the field component along the c-axis, 
H± = H cos 9. In other high-T c superconductors such as 
Bi2Sr2CaCu208, the London relation M ~ ln(r)H C 2/ H±) 
provides a good approximation for the field dependence 
of the magnetizationO Significant deviations from Lon- 
don behaviour have been noted close to T c , and TT , '- | 
preted in terms of classical (vortex) fluctuationsO'O 
Smaller deviations at lower temperatures have, been in- 
terpreted in terms of quantum fluctuations P 3 " 3 ^ As we 
now discuss, the reversible magnetization of overdoped 
Tl2Ba2Cu06+5, however, departs drastically from Lon- 
don behaviour for all temperatures and fields. 

The magnetization M(H±) can be retrieved from the 
t(9) measurements, using the relation r = MHsm9. 
This assumes that M lies close to the c-axis, which cer- 
tainly holds for vortex magnetization (but might be inval- 
idated by strong in-plane normal state paramagnetism). 
As an example, Fig. ^ shows the magnetization data on 
the T c — I5K crystal at a temperature of T ~ 1 2K, as re- 
trieved from the t{9) measurements at applied fields H 
of IT, 3T, 5T, and 10T. 

Unfortunately, the temperatures in Fig. ^ do not coin- 
cide perfectly; moreover, data points with H pointing 
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TABLE I. Compilation of the fit parameters to the experimental t(6) curves, according to Eq. (g), for all three samples and 
all fields and temperatures employed. This table provides a comprehensive representation of our experimental results. The fit 
was performed over the whole reversible region of the experimental data (except for a region of ±3° around the afe-plane where 
t(&) was too steep). The random and systematic errors in the parameters are comparable to those of individual data points, 
discussed in the Appendix 0, but the systematic error is enhanced by a factor of ~ 7 for A\ and A3, and ~ 15 for A2, due to 
the nonorthogonality of the fit functions employed in Eq. (^). For partly irreversible r(6), the random errors get worse because 
of the restricted fitting range. 



T c = 15K 


T (K) 


Ai (kN/m 2 ) 


A 2 (kN/m 2 ) 


A 3 (kN/m 2 ) 


T c = 25K 


T (K) 


A x (kN/m 2 ) 


A 2 (kN/m 2 ) 


A 3 (kN/m 2 ) 


B = 1T 


4.2 


-0.069 


0.356 


-0.028 


B = IT 


16.5 


0.32 


1.19 


-1.58 




7.3 


-0.004 


0.111 


-0.081 




19.8 


0.17 


0.51 


-0.76 




8.8 


0.031 


0.024 


-0.084 




22.9 


0.09 


0.17 


-0.22 




9.9 


0.021 


0.021 


-0.053 


B = 3T 


6.6 


5.62 


-4.22 


-15.6 




11.0 


0.009 


0.025 


-0.029 




8.3 


4.49 


-2.93 


-13.3 




12.1 


0.015 


0.009 


-0.017 




9.8 


3.67 


-2.10 


-11.4 




13.2 


0.013 


0.003 


-0.008 




11.2 


2.90 


-1.36 


-9.59 




14.3 


0.011 


0.001 


-0.004 




12.6 


2.44 


-1.13 


-8.15 


B = 3T 


4.2 


0.731 


-0.590 


-0.724 




14.1 


2.00 


-0.95 


-6.62 




7.3 


0.400 


-0.312 


-0.376 




15.6 


1.93 


-1.35 


-5.57 




8.7 


0.304 


-0.208 


-0.251 




17.0 


1.69 


-1.35 


-4.38 




9.8 


0.219 


-0.093 


-0.146 




18.5 


1.14 


-0.71 


-2.99 




10.8 


0.182 


-0.060 


-0.094 




20.0 


0.93 


-0.62 


-2.04 




11.9 


0.150 


-0.035 


-0.053 




21.5 


0.53 


-0.10 


-1.09 




13.0 


0.132 


-0.024 


-0.029 




22.9 


0.27 


0.20 


-0.48 




14.0 


0.112 


-0.008 


-0.012 




24.4 


0.22 


0.16 


-0.23 


B = 5T 


1.4 


1.46 


-0.04 


-0.31 


B = 5T 


4.2 


13.0 


-14.7 


-26.5 




2.3 


1.62 


-0.54 


-0.67 




7.0 


8.6 


-9.2 


-19.5 




4.2 


1.50 


-1.00 


-0.93 




7.7 


8.4 


-9.2 


-18.9 




7.4 


0.86 


-0.54 


-0.48 




9.3 


7.9 


-9.5 


-17.6 




8.8 


0.69 


-0.37 


-0.33 




10.7 


7.3 


-9.2 


-15.9 




9.7 


0.56 


-0.22 


-0.22 




12.1 


6.5 


-8.3 


-13.8 




10.7 


0.48 


-0.15 


-0.15 




13.5 


6.1 


-8.1 


-12.2 




11.5 


0.40 


-0.07 


-0.08 




14.9 


5.8 


-7.9 


-10.5 




12.5 


0.35 


-0.04 


-0.05 




16.3 


4.1 


-5.2 


-7.7 




13.5 


0.29 


0.01 


-0.01 




17.7 


3.2 


-3.9 


-5.7 




14.4 


0.26 


0.02 


0.00 




19.1 


2.4 


-2.7 


-4.0 


B = 10T 


1.4 


3.64 


-1.52 


-1.24 




20.6 


1.6 


-1.4 


-2.4 




2.2 


3.48 


-1.24 


-1.02 




22.0 


1.0 


-0.4 


-1.3 




4.2 


3.32 


-1.31 
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first term in Eq. (||)). This region extends as the tem- 
perature rises, and above T c , the magnetization is linear 
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FIG. 6. Diamagnetic magnetization plotted against the ap- 
plied field component H± along the c-axis, for the T c = I5K 
sample at T ~ 12K. The data were extracted from the t(9) 
measurements at applied fields H of IT, 3T, 5T, and 10T. 
For H close to the c-axis, i.e. for H ~ H±, the measurement 
uncertainty associated with a single data point diverges as 
1/ sin 0. 

close to the c-axis (i.e. for H ~ H±) have error bars di- 
verging as l/sin#. Therefore, to get a more systematic 
and comprehensive representation of our data, we em- 
ployed an averaging technique. At a given temperature, 
instead of working with individual data points, we re- 
trieved the t{6) curves from Eq. ©, where the fit coef- 
ficients A\, A2, and A3 were interpolated at the desired 
temperature from Table |. The resulting four M(H±) 
curves at the applied fields H of IT, 3T, 5T, and 10T 
were then averaged. (In order to minimize the overall 
uncertainty, a weighting factor of l/a 2 (H, H±) was as- 
signed to each curve, where <j(H,H±) is the estimated 
uncertainty of the extracted magnetization - which di- 
verges as H±_ — > H). 

This averaged magnetization is plotted against H± and 
logH± in Fig. [?] for a variety of temperatures below T c 
for all three crystals. The dashed lines in all the pan- 
els represent estimates for the equilibrium magnetiza- 
tion in fields below the irreversibility field. These esti- 
mates were obtained by extrapolating the fit in Eq. (j^) 
to angles higher than 8[ rr and therefore to applied fields 
H± < H ilr . The reader should ignore the small "wiggles" 
in the lines: these are artefacts, caused by slightly differ- 
ent values of M obtained from runs in the experimentally 
applied fields of IT, 3T, 5T, and 10T. 

Consider the right panels in Fig. 0, showing the varia- 
tion of the magnetization with log H± . For the T c = 85K 
crystal, there are some departures from linearity, but the 
form of the data is qualitatively similar tq_ that seen 
in other cuprates such as Bi2Sr2CaCu2O8.ES For the 
T c = 15K and 25K crystals, however, M(\ogH±) departs 
strongly from linearity and displays an upturn which - 
at least at lower temperatures - occurs in the vicinity of 
the irreversibility field. As can be seen in the left panels 
of Fig. |7|, there then exists an extensive field range over 
which M ~ H± (corresponding to contributions from the 



over the whole field range, as discussed in section |IVC 

A (roughly) linear magnetization, M ~ H±, has been 
observed previously in other high-T c materials in the im- 
mediate vicinity of T c , and ascribed to classical thermal 
fluctuations.L-3 Its appearance at very much lower tem- 
peratures is a surprising result and raises the possibility 
that some sort of artefact might be affecting the mea- 
surements. A linear relation between M and H is char- 
acteristic of both the Pauli susceptibility in a normal 
metal and a linear Curie spin susceptibility. All high- 
Tc materials exhibit such a spin susceptibility which is 
anisotropic due to crystal field interactions,^ leading to 
a torque that varies as H 2 sin 28. Such a (normal state) 
torque mimics the first term in Eq. (^|), suggesting that 
the linear magnetization may have nothing to do with 
superconductivity and may simply reflect the presence 
of a normal-state susceptibility anisotropy with a strong 
temperature dependence. 

Although possible in principle, there are two reasons 
why such a scenario is unlikely: firstly, the sign of the 
effect is implausible. Fig. || shows that the sin 28 com- 
ponent changes sign at higher temperatures, for both 
the T c = 15K and 25K crystals. At higher temperatures, 
the magnitude of the associated susceptibility anisotropy 
is ~ 10~5_.(SI units), a number that is typical for the 
cuprates ;E3 the sign of this term is such that the maxi- 
mum susceptibility lies along the c-axis, as observed in es- 
sentially all other high-T c materials. However, as shown 
in Fig. H, the sign is reversed below T ~ 2.5T C , for both 
crystals. It is very difficult to understand why a nor- 
mal state effect of this type should flip sign, and why it 
should do so at a temperature linked to the superconduct- 
ing transition. In addition, at a given temperature, the 
normal state susceptibility in Ti2Ba2Cu06+,5 is known ta 
be larger for crystals with lower critical temperatures. 
The linear susceptibility deduced from our torque mea- 
surements is shown in Fig. |^; at a given temperature, it 
is smaller for the lower T c crystal. Although a normal 
state origin for the linear magnetization cannot be com- 
pletely dismissed, the balance of evidence suggests that 
it is likely to be a consequence of superconducting order. 

It should be mentioned that the diamagnetic upturn 
in our susceptibility anisotropy data in Fig. ^| exhibits 
qualitative similarity to the "diamagnetic Curie- Weiss, 
law" observed in the normal state of YBa2Cu307l23€£l 
and, most notably, in La2- £C Sr ;! ;Cu204,cEI the origin of 
which remains unknown. Fitting our A%(T) data to a 
Curie- Weiss term with a constant background, we obtain 
025K — 15K and, rather intriguingly, 815K — as char- 
acteristic temperatures for the two samples, where the 
subscript refers to the critical temperature. 

It is also interesting to note that, independent of 
the origin of the sin 28 contribution in overdoped 
Tl2Ba2Cu06+5, subtracting it from the data does not 
lead to a standard superconducting torque. Although its 
angular dependence is qualitatively described by Eq. (j|), 
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FIG. 7. Reversible diamagnetic magnetization (solid lines) and its extrapolation into the irreversible region (dashed), plotted 
against the applied field component H± along the c-axis (left panels) and its logarithm (right panels) , at different temperatures 
for all three samples. The data were interpolated from the fit parameters in Table |, assuming that the magnetization is always 
parallel to the c-axis. The estimated error is around ±2A/m for the T c — 15K sample, ±15A/m for T c — 25K, and ±50A/m 
for T c — 85K. The "wiggles" in the lines are artefacts, caused by slightly different values of M obtained from runs in the 
experimentally applied fields of IT, 3T, 5T, and 10T. 
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attempts to obtain the apparent upper critical field i]H C 2 
gave unphysical temperature dependences. This is not 
the case for other high-T c superconductors, where physi- 
cally reasonable values and temperature-dependences of 
H c2 can be extracted from torque data.cJH 

If all the data shown in Fig. [?] are indeed directly linked 
to superconducting order, they give a new perspective on 
one of the problems that motivated our study, that of the 
upper critical field in strongly overdoped Tl2Ba2Cu06+<5. 
Recently, an attempt has been made to deduce the tem- 
perature dependence of H C 2 from a Raman scattering ejsu 
periment on a Tl2Ba2Cu06+<5 crystal with T c ~ 25K.E3 
In our T c = 25K sample, our data suggest that diamag- 
netism persists to much higher fields than those identified 
as H C 2 in that study. 

Our data also clearly imply the existence of dia- 
magnetism - and thus remnant superconducting or- 
der - at fields and temperatures for which transport 
measurements^ suggest that the normal state has been 
fully restored. This conclusion was previously reached 
by Carrington_et al. on the basis of heat capacity 
measurements,!^ but the field dependence of the heat ca- 
pacity anomaly in low-T c T^E^CuOg+s appeared to be 
qualitatively similar to that seen in other cuprates with 
higher transition temperatures. As noted in Ref. |2^, this 
is difficult to reconcile with striking differences in the 
transport behaviour between optimally doped and over- 
doped Tl2Ba2Cu06+5. The new result from our torque 
data is that although remnant superconducting order ex- 
ists above the irreversibility line for all doping levels, de- 
partures from London behaviour clearly emerge as T c is 
reduced from 85K to 15K. 

In summary, the magnetic behaviour of low-T c crys- 
tals belonging to the Tl2Ba2Cu06+<5 system has been a 
puzzle for a number of years. Transport measurements 
have suggested that a vortex liquid does not appear and 
that the irreversibility and upper critical field bound- 
aries coincide. On the other hand, specific heat measure- 
ments have shown that superconducting order persists 
well above the irreversibility boundary. In this work, 
we have studied the superconducting magnetization and 
demonstrated that significant departures from London 
liquid behaviour occur for all three samples studied, i.e. 
across the entire overdoped region. The departures in- 
crease in severity as the doping level is increased. In 
particular, for strongly overdoped (T c = 15K) material, 
a London-like vortex liquid does not appear above the 
irreversibility boundary. The remnant superconducting 
order is characterized by a linear diamagnetic response, 
a response that persists well above T c and also up to the 
highest field employed here (10 Tesla). We hope that this 
novel result will encourage further effort to develop our 
understanding of the intriguing behaviour of overdoped 
Tl 2 Ba 2 Cu0 6+5 . 
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APPENDIX A: CALIBRATION 

With a force constant of 2.5 N/m, the angular displace- 
ment of the lever for all the reversible torques encoun- 
tered in this study is estimated to be significantly less 
than 1°. Under these circumstances, it can be assumedo 
that the change in lever resistance, AR, is proportional 
to the torque r experienced by the sample, or, 



t = aAR 



(Al) 



where a is a calibration constant. 

Let p and V be the density and volume of the attached 
crystal, respectively. To calibrate the cantilever, the bal- 
anced bridge configuration described in the text was ro- 
tated in zero magnetic field where the lever only experi- 
ences a gravitational torque given by 



t„ = p Vg I sin 9 



(A2) 



where 9 is the angle to the vertical, I is the length of the 
lever, and g is the gravitational acceleration of the Earth. 
The resistive path is assumed to be close to the region of 
maximum stress at the lever base. The amplitude of the 
associated resistance change, AR g , was measured and the 
sensitivity obtained from 



pvgi 

AR a 



(A3) 



The value of a obtained in this way was 
8 x 10 _11 Nm/f2. No systematic study of the tempera- 
ture dependencej-Oif the sensitivity was performed. How- 
ever, Yuan et al.L3 found the piezolever sensitivity to be 
constant to a factor of 2 over the temperature range in- 
volved in this work, a result that was confirmed for one 
lever used in this study. The calibration constant in the 
present work is therefore subject to an uncertainty of or- 
der ±50%. 



APPENDIX B: ERROR ESTIMATE 

The main source of random measurement uncertainty 
stemmed from temperature variations. Although these 







variations were negligible when the lever was immersed 
in liquid helium, they were about ±0.03K when it was 
surrounded by exchange gas. These variations coupled 
into the system via the temperature dependent lever re- 
sistance R(T). The compensation arrangement could not 
eliminate this effect, but at least suppress it to 5% or 
less. Since dR/dT itself varied, the resulting torque un- 
certainty depended on temperature and rose from to a 
minimum of about 2 x 10 -13 Nm at 5K to a maximum 
of about 10" 12 Nm at 30K. 

The small contribution from the gravitational signal 
(Appendix |A|) was measured and subtracted from the 
data. Cantilever magnetoresistance was less easy to deal 
with and constituted the most serious systematic uncer- 
tainty for the measurements reported here. Initially vary- 
ing as B 2 , the (temperature-dependent) change in lever 
resistance with field becomes almost linear above 0.5T 
and reaches a value as high as ~ 100 fl in a field of 10T. 
Again, the compensation lever could trim down the effec- 
tive change to a few Ohms, but this was still comparable 
to the total signal for the (T c — 15K) crystal. 

The strategy for reducing this uncertainty was two- 
fold. In the first place, although the magnetoresistance 
is large, its anisotropy at temperatures below 40K is 
relatively small, reducing the apparent magnetoresistive 
torque to 6 x 10~ n N/m 2 T when the levers were rotated 
in the field. By confining measurements to the angu- 
lar variation of the torque, the magnetoresistance error 
was therefore reduced to a manageable level. (This was 
the reason why the irreversibility field was obtained from 
the angular characteristic, rather than a traditional field 
sweep.) 

It was further reduced by noting that, in our geomet- 
rical configuration, the symmetry of the vortex torque r v 
and the apparent magnetoresistive torque r m around the 
ab-plane were odd and even, respectively: 

T V (7T - 0) = -T v (0) (Bl) 

r m (7r - 0) = r m (6) 

Given an observed angular dependence r raw (0), over 
two angular quadrants, the vortex torque t v {9) was there- 
fore extracted according to: 

ry(0) = \ - r raw (^ - 6)) (B2) 

This is the quantity actually plotted in Figs. § and ^. 
This procedure cancelled the symmetric magnetoresis- 
tive component. Misalignments of the levers (~1°) led 
to a residual torque uncertainty of ~ 10~ 13 Nm/T (be- 
low 40K; ~ 10~ 12 Nm/T at higher temperatures). Note 
that, in a 10T field, this is about 10 2 times larger than the. 
intrinsic uncertainty quoted for the piezolever method.Ej 
The degree of cancellation could probably be improved, 
but it was adequate for the present work. In a field of 
10T, the remaining torque uncertainty corresponds to a 
moment uncertainty of ~ 10 _13 Am 2 , about 10 2 times 



better than can be achieved by commercial SQUID meth- 
ods. 

To compare the above numbers with the torque density 
results, they have to be divided by the respective crystal 
volumes Vi 5K = 2.7 x 10~ 13 m 3 , U 25 k = 1.3 X 10" 13 m 3 , 
and V85K = 0.65 x 10 _13 m 3 , where the subscript refers 
to the critical temperature. 
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